Comparison of two non - primitive methods for path integral simulations : Higher - order corrections vs . an effective propagator approach

Two methods are compared that are used in path integral simulations. Both methods aim to achieve faster convergence to the quantum limit than the so-called primitive algorithm (PA). One method, originally proposed by Taka-hashi and Imada, is based on a higher-order approximation (HOA) of the quantum mechanical density operator. The other method is based upon an effective propagator (EPr). This propagator is constructed such that it produces correctly one and two-particle imaginary time correlation functions in the limit of small densities even for finite Trotter numbers P. We discuss the conceptual differences between both methods and compare the convergence rate of both approaches. While the HOA method converges faster than the EPr approach, EPr gives surprisingly good estimates of thermal quantities already for P = 1. Despite a significant improvement with respect to PA, neither HOA nor EPr overcomes the need to increase P linearly with inverse temperature. We also derive the proper estimator for radial distribution functions for HOA based path integral simulations.